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Energy conservation via the pathway of dissimilatory sulfate reduction is present in a 
diverse group of prokaryotes, but is most comprehensively studied in Deltaproteobacteria. 
In this study, whole-genome microarray analyses were used to provide a model of 
the energy metabolism of the sulfate-reducing archaeon Archaeoglobus fulgidus, based 
on comparative analysis of litoautotrophic growth with H2/CO2 and thiosulfate, and 
heterotrophic growth on lactate with sulfate or thiosulfate. Only 72 genes were expressed 
differentially between the cultures utilizing sulfate or thiosulfate, whereas 269 genes 
were affected by a shift in energy source. We identified co-located gene cluster encoding 
putative lactate dehydrogenases (LDHs; lldD, did, lldEFG), also present in sulfate-reducing 
bacteria. These enzymes may take part in energy conservation in A. fulgidus by specifically 
linking lactate oxidation with APS reduction via the Qmo complex. High transcriptional 
levels of Fqo confirm an important role of F420H2, as well as a menaquinone-mediated 
electron transport chain, during heterotrophic growth. A putative periplasmic thiosulfate 
reductase was identified by specific up-regulation. Also, putative genes for transport of 
sulfate and sulfite are discussed. We present a model for hydrogen metabolism, based 
on the probable bifurcation reaction of the Mvh:Hdl hydrogenase, which may inhibit the 
utilization of Fd^Q^fo^ energy conservation. Energy conservation is probably facilitated via 
menaquinone to multiple membrane-bound heterodisulfide reductase (Hdr) complexes 
and the DsrC protein — linking periplasmic hydrogenase (Vht) to the cytoplasmic reduction 
of sulfite. The ambiguous roles of genes corresponding to fatty acid metabolism induced 
during growth with H2 are discussed. Putative co-assimilation of organic acids is favored 
over a homologous secondary carbon fixation pathway, although both mechanisms may 
contribute to conserve the amount of Fd^ed needed during autotrophic growth with H2. 

Keywords: Archaeoglobus fulgidus, hydrogenase, heterodisulfide reductase, dissimilatory sulfate reduction, 
lactate dehydrogenase 



INTRODUCTION 

The sulfate-reducing prokaryotes (SRP) have played a central 
role in cycling of carbon and sulfur in anoxic environments 
throughout long periods of Earth's geological history. Despite 
early characterization of the cytoplasmic pathway of dissimilatory 
sulfate reduction (Peck, 1962) it is only in recent years that the 
mechanisms facilitating energy conservation in SRP have been 
more comprehensively characterized (Pereira et al., 2011; Grein 
et al.,2013). 

The genus Archaeoglobus comprises of archaeal, 
(hyper) thermophilic, dissimilatory sulfate reducers (Stetter et al., 
1987; Stetter, 1988) and is phylogenetically associated with the 
lineages of Methanosarcinales, Methanomicrobiales, and uncul- 
tured ANME-1 (Brochier-Armanet et al, 2008; Guy and Ettema, 
2011). The type species A. fulgidus VC16 is a chemolithoau- 
totroph that utilizes H2 or formate as electron donors for 
autotrophic growth. In addition, A. fulgidus grows carboxy- 
dotrophically on CO/CO2 and as a chemoorganoheterotroph 
utilizing a wide range of substrates including fatty acids, alkenes. 



complex peptides, and specific amino acids (Stetter et al., 1987; 
Stetter, 1988; Hartzell and Reed, 2006; Henstra et al, 2007; 
Khelifi et al, 2010; Parthasarathy et al, 2013). For the complete 
oxidation of organic substrates to CO2, A. fulgidus uses a modi- 
fied acetyl- Co A pathway with similar enzymes and cofactors as in 
the methanogens (Moller-Zinkhan et al., 1989; Moller-Zinkhan 
and Thauer, 1990; Vorholt et al, 1995; Estelmann et al, 2011). 
Reduction of sulfate (S04^~) to sulfide (S^~) in A. fulgidus 
proceeds via the highly conserved dissimilatory sulfate reduction 
pathway of the SRP (Peck, 1962; Klenk et al., 1997; Pereira 
et al., 2011). This was probably acquired by Archaeoglobales via 
multiple lateral gene transfer events from an early ancestor of 
clostridial SRP (Klein et al., 2001; Zverlov et al, 2005; Meyer and 
Kuever, 2007). 

The energy conservation mechanisms in A. fulgidus are incom- 
pletely understood. During growth on lactate, the reduced coen- 
zyme F420 (F420H2) is generated from the oxidative acetyl-CoA 
pathway. The presence of both menaquinone and a homolog 
of the respiratory NAD(P)H:quinone oxidoreductase complex. 
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the F42oH2:quinone oxido reductase complex (Fqo), suggest that 
electrons from F420H2 are transferred to the membrane-bound 
respiratory chain by the Fqo complex. Fqo probably couples the 
reduction of menaquinone and proton translocation. (Tindall 
et al., 1989; Kunow et al., 1993; Baumer et al., 2000; Brliggemann 
et al., 2000) A D-lactate dehydrogenase is confirmed to be present 
(Reed and Hartzell, 1999), but it is unclear how this membrane 
associated enzyme facilitates energy conservation, as it is shown 
to interact with a NADH oxidase (Pagala et al., 2002). Also, the 
CO factor NAD(P)H plays a negligible role in energy conserva- 
tion (Noll and Barber, 1988; Kunow et al, 1993; Warkentin et al, 
2001). 

There is also a possible alternative energy conservation path- 
way in A. fulgidus. In D. vulgaris, cytochrome c mediated "hydro- 
gen cycling" is suggested as an energy conservation mechanism 
during growth with lactate (Odom and Peck, 1981; Keller and 
Wall, 2011). In this reaction, formation of hydrogen is a result 
of cytoplasmic oxidation of lactate. The subsequent diffusion 
and periplasmic oxidation of hydrogen contributes to the for- 
mation of a proton gradient. In Methanosarcina barkeriy the 
Vht/Vhx dehydrogenase also facilitates a hydrogen cycling mecha- 
nism under heterotrophic growth conditions, and sustains growth 
when Fpo (Fqo) is absent in deletion mutants (Kulkarni et al, 
2009). The presence of a cytoplasmic as well as a periplasmic 
hydrogenase in A. fulgidus (Mander et al., 2004) potentially fulfills 
requirements for a "hydrogen- cycling" mechanism. 

Two co-located heterodisulfide reductase (Hdr) -associated 
hydrogenases are present in the genome of A. fulgidus, which 
are homologous to those involved in energy conservation in 
the methanogens (Mander et al, 2004). These are the solu- 
ble [NiFe]hydrogenase/heterodisulfide-like (MvhABC/HdlABC) 
complex and the membrane-bound uptake hydrogenase, "F420- 
non-reducing hydrogenase" (Vho/Vht). Reduced ferredoxin 
(Fdred) is essential for fixation of CO2 through the acetyl- 
CoA pathway. In methanogens, the Mvh:Hdl complex homolog, 
Mvh/Hdr, couples the exergonic reduction of the heterodisulfide, 
CoM-S-S-CoB, with endergonic reduction of ferredoxin with H2, 
by a flavine-based bifurcation mechanism (Kaster et al., 2011). 
The periplasmic Vht hydrogenase reduces the quinone-like cofac- 
tor methanophenazine coupled to the membrane-bound HdrDE, 
facilitating energy conservation during growth on H2 (Ide et al, 
1999; Thauer et al, 2010). 

Despite the absence of genes and cofactors for terminal 
methanogenesis (Stetter et al., 1987; Klenk et al., 1997), sev- 
eral factors suggest that thiol/disulfide conversions catalyzed by 
Hdr are involved in electron transfer and energy conservation 
in A. fulgidus, as has been proposed for methanogens and more 
recently for SRP (Mander et al, 2002, 2004; Pereira et al, 2011; 
Grein et al, 2013). All known SRP, including A. fulgidus, encode 
HdrA and HdrDE related genes, which almost ubiquitously form 
membrane-bound redox complexes (Pereira et al., 2011; Grein 
et al., 2013). These complexes may facilitate energy conserva- 
tion during different steps of sequential dissimilatory sulfate 
reduction. The quinone-interacting membrane-bound oxidore- 
ductase (QmoABC) complex probably links the electron transfer 
chain to the first reductive step of sulfate reduction catalyzed 
by adenosine- 5^ -phosphosulfate (APS) reductase (AprAB) (Pires 



et al., 2003; Zane et al., 2010; Grein et al., 2013). In Desulfovibrio 
it has recently been proposed that the Qmo subunit homologous 
to the bifurcating HdrA, QmoB, may facilitate a "confurcation" 
mechanism (Ramos et al., 2012). The "confurcating" Qmo com- 
plex may catalyze energy conservation by proton translocation 
via an endergonic periplasmic menaquinol oxidation, driven by 
an exergonic cytoplasmic oxidation reaction coupled to terminal 
reduction of APS. The second complex, DsrMK, is a homolog of 
HdrDE, and is ubiquitous amongst SRP (Pereira et al., 201 1). This 
complex probably facilitates energy conservation and is linked by 
electron transfer via disulfide/thiol redox reactions, to the termi- 
nal step of sulfite reduction by bisulfite reductase/sulfite reduc- 
tase (DsrAB) (Mander et al, 2002; Pires et al, 2006). Similarly 
to the HdrDE of methanogens, the DsrMK complex probably 
couples periplasmic oxidation of reduced menaqinone (instead 
of reduced methanopenazine) to cytoplasmic cysteine disulfide 
(Cys-S-S-Cys) reduction, in the enzyme DsrC (in stead of a CoM- 
S-S-CoB) (Mander et al, 2005). Unusually, dsrMK is encoded 
by multiple homologs in A. fulgidus, corresponding to multiple 
DsrMK and a DsrMK(JOP) complex, which differ in domain 
composition and among lineages of SRP (Klenk et al., 1997; 
Pereira et al, 201 1). The dsrC gene is ubiquitously present in SRP, 
and DsrC is the probable link between heterodisulfide reductase 
(DsrK) and DsrAB (Oliveira et al, 2008; Pereira et al, 201 1; Grein 
et al., 2013). However, it should be noted that although it is likely 
that the DsrMK(JOP) complexes may facilitate proton transloca- 
tion by MQH2 oxidase:DsrC reductase, it is questioned whether 
this reaction is thermodynamically favorable (Thauer et al., 2007; 
Grein et al, 2013). 

The role of reduced ferredoxin (Fdred) in energy conservation 
in SRP remains unclear, as it has been proposed as an electron 
donor for both APS and sulfite reduction (Oliveira et al., 2008, 
2011; Ramos et al., 2012). In A. fulgidus, this offers a potential 
coupling between ferredoxin and electron transport phosphory- 
lation, but also represents a significant bioenergetic challenge, as 
fixation of CO2 through the acetyl-CoA pathway requires Fdred- 
Interestingly, while chemoorganotrophic and carboxydotrophic 
growth are coupled to sulfate reduction in A. fulgidus, only thio- 
sulfate or sulfite are utilized with H2 as energy source (Stetter 
et al., 1987; Steinsbu et al., 2010). This may potentially be coupled 
to the role of Fdred in energy and carbon metabolism. 

To provide a deeper insight into electron transport and energy 
conservation mechanisms in A. fulgidus, we used whole genome 
microarrays to identify redox complexes expressed under differ- 
ent growth conditions. Previously, only the heat shock response in 
A. fulgidus has been characterized by global transcriptional pro- 
filing (Rohlin et al., 2005). We examined heterotrophic growth 
with lactate and litoautotrophic growth with H2, as well as the 
differential use of the electron acceptors thiosulfate and sulfate. 
The results form an overall energy conservation model where 
the Fqo and membrane-bound electron transport, facilitated by 
menaquinone, Qmo and multiple DsrMK, are central to energy 
conservation during growth with lactate. During growth with 
hydrogen, our model suggests that Fdred ? generated by Mvh:Hdl, 
is utilized primarily for carbon assimilation and probably does 
not contribute to energy conservation. From the data and com- 
parative genomics it seems likely that the inability of A. fulgidus 
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to grow with sulfate when hydrogen is an energy source is caused 
by transcriptional regulation of the gene for pyrophosphatase, 
resulting in the blocking of APS formation. Overall, the results 
point to a key role in energy conservation for electron trans- 
fer from hydrogen to thiosulfate, facilitated by thiol/disulfide 
conversions catalyzed by membrane-bound DsrMK in A. fulgidus. 

METHODS 

Archaeoglobus fulgidus strain VC16 (DSMZ 4302) obtained from 
the Deutsche Sammlung von Mikroorganismen und Zellkulturen 
(Braunschweig, Germany) was cultivated in anoxic, carbonate 
buffered medium (10 ml medium in 26 ml serum vials) under an 
atmosphere of N2:C02 80:20 (1 atm), at pH 6.8. The composition 
of the media was as follows: 0.32 g/1 KCl, 1.0 g/1 MgCl2«6H20, 
0.25g/l NH4CI, 0.14g/l CaCl2«2H20, 0.11 g/1 K2HP04-3H20, 
0.2 g/1 KH2PO4, 18.0 g/1 NaCl, and 0.3 g/1 yeast extract. Minor 
constituents were; 0.015 g/1 Titriplex I (Nitriloaceticacid), 
0.005 g/1 MnS04»2H20, 0.001 g/1 CoCl2«6H20, 0.001 g/1 
ZnS04«7H20, 0.0001 g/1 CuS04«5H20, 0.0001 g/1 H3BO3, 
0.0001 g/1 Na2Mo04«2H20, 0.002 g/1 NiS04«6H20, 0.039 mM 
(NH4)2Fe(S04)2«6H20, and 0.5 ml/1 0.2% Resazurin. After 
autoclaving, sterile anoxic solutions were added to the medium 
to a final concentration of 30 mM NaHCOsand 0.25 mM Na2S. 
When thiosulfate was used as an electron- acceptor, sulfate; 
2.2 g/1 Na2S04 and 3.7 g/1 MgS04«7H20, was exchanged with 
3.7 g/1 MgCl26H20 and 7.45 g/1 Na2S203«5H20 (thiosulfate). All 
media contain 0.18 mM S04^~ attributed to the composition of 
minor constituents. Filter- sterilized (0.2 |xM), anoxic, thiosulfate 
solution was added after autoclaving. During heterotrophic 
growth, 35 mM sodium-D,L-lactate (50/50) was added to the 
medium whereas 250 kPa H2:C02 (80:20 ratio) was used during 
litoautotrophic growth. 

Cultivation was performed at 80° C and the tubes were incu- 
bated at an approximate angle of 6° ensuring a high surface to 
volume ratio. The turbidity of samples (absorbance at 600 nm) 
was used for monitoring of cultures; linearity against direct cell 
counts for all growth conditions was confirmed throughout the 
absorbance range, (using a Thoma- chamber; depth 0.02 mm). All 
cultures were harvested at a pre-determined absorbance, on the 
basis of growth experiments. 

Cultures were flash cooled (20 s) to approximately 0°C in 
a — 80°C, 70% ethanol: water slurry, and harvested by centrifuga- 
tion in 15 ml falcon tubes at 3000 g for 15 min at 0°C. The pellet 
was immediately re-suspended in 100 |xl RLT buffer of the RNeasy 
kit (Qiagen) and stored at — 80°C for a maximum of 1 week before 
total RNA extraction. 

PREPARATION OF RNA 

Upon RNA extraction, samples were thawed to 37°C and placed 
on ice. Samples from equivalent growth conditions were pooled 
in numbers sufficient to obtain the required yield (1 |xg). The 
final volume was adjusted to 600 |xl with buffer (RLT, RNeasy) 
before proceeding. The RNeasy mini kit (Qiagen) was used 
for the total RNA extraction and with an additional DNase I 
(Qiagen) step. RNA was concentrated using the RNeasyMinElute 
kit (Qiagen) to achieve the required concentration of 1 |xg/ml 
total RNA for cDNA synthesis. Total RNA concentration was 



determined photometrically (Cary 300 UV-Vis, Varian) using a 
TrayCell cuvette with a 0.2 mm cap (Hellma, Germany), and 
RNA quality was evaluated using a RNA 6000 Nano kit with a 
Bioanalyzer 2100 instrument (Agilent). 

Each sample prepared for hybridization was the result of pool- 
ing in order to obtain sufficient material for analysis. Samples 
consist of 5 individual randomly selected tubes for T-H2/CO2 
samples, and 2 individual tubes for all S-L and T-L samples. A 
total of 27 microarray hybridizations were performed. These cor- 
respond to the following growth conditions (outline of design; 
Figure lA): S-L, 10 hybridizations (mid-log 6, late log 4); T-L 
6 hybridizations (mid-log); T-H2/CO2 11 hybridizations (mid- 
log 7, late log 4). The arrays utilized were the commercially 
available Roche Nimblegen 080626 Aful DSM4304 design, uti- 
lizing the 4-plex array design (4 x 72K format) where 2392 
open reading frames (ORF's) are assayed. RNA was prepared 
according to the guidelines provided for the Nimble Chip arrays 
(Nimblegen systems, 2007); double stranded cDNA was syn- 
thesized by Superscript Double-Stranded cDNA synthesis kit 
(Invitrogen), with Random Hexamer Primer (Roche Applied 
Science). Subsequent steps were performed as recommended by 
the array supplier. 

The analysis of the obtained image, and quality control 
was performed by the NimbleScan software version 2.5, values 
obtained from robust multiarray average (RMA) normalization 
(Irizarry et al, 2003) was utilized for further analysis. The data 
were deposited in the ArrayExpress database (https://www.ebi.ac. 
uk/arrayexpress/) under the accession code: E-MTAB-2294. 

RMA normalized data were further analyzed by the J-Express 
software 2012 (http://jexpress.bioinfo.no/site/), and quantile nor- 
malization was performed on all samples (Bolstad et al, 2003). 
Due to the quantile normalization procedure, the mean inten- 
sity of all assayed genes is practically constant. Therefore, the 
mean signal intensity was arbitrarily set as 1.0 and values are 
reported relative to this level in order to convey the level of signal 
abundance. 

Analysis of variance (ANOVA) is the principle method of sta- 
tistical evaluation in this work. In order to perform ANOVA, a 
similar sample size is required. The 27 samples correspond to five 
different growth conditions and were of unequal size. To increase 
robustness of analysis, the minimum residual sum of squares per 
gene from the total data was selected for ANOVA. Therefore, 4 
values representing each of the corresponding growth conditions 
were analyzed; S-L, late-log S-L, T-L and T-H2/CO2, and late-log 
T-H2/CO2. We report differential expression as significant when 
the ANOVA returns ap-value of less than 0.00001 {p < 0.00001, 
critical F > 37.71). A significant shift between two groups cor- 
responding to 1.5 fold or larger was generally considered as a 
cut-of for major differential expression. Correspondence analy- 
sis was used for clustering of individual samples (Fellenberg et al, 
2001) and as a control of the ANOVA (Figure IB). 

Functional annotation was performed using the latest version 
of archaeal clusters of orthologous genes (arCOG) (Wolf et al, 
2012). Enrichment analysis was performed using analysis between 
selected groups and entire dataset using the Chi- squared test (p < 
0.05 for groups larger than 5 genes). Association of individual 
genes to KEGG pathways were retrieved from the KEGG database 
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FIGURE 1 I Growth curves and illustration of differentially regulated 
genes. (A) Representative growth curves of cultures grown with lactate and 
sulfate (S-L), lactate and thiosulfate (T-L), hydrogen/C02 and thiosulfate 
(T-H2). Boxes indicate absorbance values when sannpling was perfornned. The 
nunnbers indicate annount of replicate nnicroarray hybridizations for each 
growth condition. (B) Correspondence analysis (CA) provides an overview of 
variance and differentially regulated genes determined by ANOVA. Genes 



(circles) and individual hybridizations (squares) are clustered with respect to 
each other by CA, and are displayed along the two first principle connponents 
(Total variance retained 41 %; Principle Connponent Axis — gray lines; PCA1 : 
27.2%, PCA2: 13.8%). Colored lines indicate nnedians corresponding to each 
growth condition. Significantly differentially expressed genes (ANOVA) are 
colored according to growth condition; black lines denote a fold change 
greater than 1.5. 



(http://www.kegg.jp/). In order to evaluate the vaUdity of signal 
intensity vs. functional genes, enrichment of genes correspond- 
ing to KEGG pathways were evaluated by a Kolmogorov-Smirnov 
statistic equivalent; Gene Set Enrichment Analysis (Subramanian 
et al, 2005). Enrichment of KEGG pathway associated genes was 
performed on a list sorted by minimal recorded signal- intensity 
per gene from any sample. 

Homology searches were conducted using either BLASTp or 
PSI-BLAST using default settings (http://blast.ncbi.nlm.nih.gov/ 
Blast. cgi), while conserved domains were identified using the 
Conserved Domains Database (CDD) database (http://www. 
ncbi . nlm . nih . go v/S tr uctur e/cdd/cdd . shtml) ( Marchler- Bauer 
et al, 201 1). Shared synteny of gene clusters were identified using 
the STRING database (http://string-db.org/), in combination 
with the de-novo synteny explorers Absynte and Syntax (http:// 
archaea.u-psud.fr/archaeasoftwarepage.html) (Despalins et al., 
2011; Oberto, 2013); which also provide the graphical templates 
for Figures 3, 4. 

RESULTS 

GROWTH AND GENERAL TRANSCRIPTIONAL SHIFT 

Doubling time was more than halved in A. fulgidus when thio- 
sulfate (T-L) was added as terminal electron acceptor instead of 
sulfate (S-L) during growth with lactate (Figure lA). The specific 
growth rate (|x) evaluated by 7 replicate cultures, was significantly 
higher for T-L (|x: 0.28 =b 0.07 h"^; doubling time 1.1 =b 0.4 h), 
than both S-L (0.12 zb 0.004h-i; 2.4 zb 0.1 h) and T-H2/CO2 
(0.13 zb 0.03 h"^; 2.3 zb 0.9 h). Turbidity was estimated to increase 
by 2.6 absorbance units per cell (A6oonm) in cultures grown with 
lactate (S-L, T-L) compared to the T-H2/CO2 cultures. Hence, 
growth rate increased during growth with T-L, vs. indistinguish- 
able rates between S-L and T-H2/CO2 cultures. The increase in 



growth yield inferred from absorbance in lactate -grown cultures 
may be partially explained by an increase in cell size during 
growth with lactate. 

The 27 transcriptional profiles of A. fulgidus cells cultivated 
with S-L, T-L, T-H2/CO2, and late log-phase cells cultivated 
with S-L and T-H2/CO2 were compared by microarray analy- 
sis (Figure lA). Correspondence analysis (Fellenberg et al., 2001) 
revealed that individual samples clustered together, with samples 
from similar growth conditions being distinct from other assayed 
conditions (Figure IB). Between any of the assayed conditions 
a total of 1268 genes were differentially expressed (ANOVA p < 
0.00001); of these 514 genes were differentially expressed over 1.5 
fold (53%; and 21% of total assayed genes). 

The analysis identified 692 differentially regulated genes (29% 
of assayed genes. Figure IB) corresponding to either electron 
donor/carbon source (S-L/T-L vs. T-H2/CO2), electron accep- 
tor (S-L vs. T-L, T-H2/CO2) or growth phase (log; S-L, T-L, 
T-H2/CO2 vs. late log; S-L, T-H2/CO2). Of these, 369 genes 
were over 1.5 fold differentially expressed (15% of assayed genes. 
Figure IB). 

The differentially regulated genes were as follows (Figure IB): 
hydrogen vs. lactate 419 genes; 264 over 1.5 fold (17.5%; 11%: 
167/97 up/down), thiosulfate vs. sulfate 121 genes; 68 over 1.5 
fold (5.1%; 2.8%: 21/47), late vs. mid-log; 192 genes; 61 over 
1.5 fold (8%; 2.6%: 33/28). Hence, a shift in energy metabolism 
and carbon source introduced the largest number of differen- 
tial regulation. Whereas, a shift in electron donor contributes 
less to the total magnitude of transcriptional regulation, at com- 
parable levels to differential expression related to growth phase. 
Differential expression corresponding to either late log S-L, late 
log T-H2/CO2, or T-L were not considered further (21% of 
assayed genes and 6% regulated above 1.5 fold). 
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FIGURE 2 I Distribution and enrichment of COG categories that are 
differentially regulated corresponding to specific factors (Figure 1B). 

Bar plots display expression as percentage of each COG category. Bars 
highlighted with specific colors (as in Figure IB) correspond to 



significantly enriched COG categories of genes induced by this factor 
(Chi-squared test; p > 0.05, nnore than 4 genes). The nunnbers denoted 
in brackets are the nunnber of differentially regulated genes in each 
group. 



DIFFERENTIALLY EXPRESSED GENES WITH HYDROGEN AND CQ2 vs. 
LACTATE 

Of the 206 genes up-regulated in cultures grown with T-H2/CO2 
vs. S-L and T-L, 92 were affiliated with the COG's corresponding 
to metabolic processes: where 3 1 genes corresponding to energy 
production (C) and 34 to lipid transport and metabolism (I) were 
significantly enriched (Figures IB, 2). There was also a major dif- 
ferential expression involving 5 genes related to cellular motility 
(N) and signal transduction (T). Specific genes in the energy pro- 
duction category included the 10 co-located hydrogenase genes 
of two distinct complexes (AF1372-AF1377, AF1379-AF1381, 
Table ID; Mander et al., 2004), and a membrane-bound Hdr 
(AF0755, Table IC). 

Surprisingly, the genes corresponding to D -lactate dehydro- 
genase (AF0394) and a putative d-lactate dehydrogenase gene 
(AF0868) (Reed and Hartzell, 1999; Pagala et al, 2002) were 
induced during growth with T-H2/CO2. The gene AF0394 
was up-regulated by a minor fold (<1.5) and co-regulated 
(Pearson s r-value; 0.89) with the associated NADH oxidase gene 
(AF0395, >1.5 fold) (Table IE, Figure 3). 

Other highly induced genes classified under energy pro- 
duction belong to the TCA cycle; succinate dehydroge- 
nase (AF0682-AF0683) and malate oxidoreductase (AF1727) 
which probably has a non- energetic, assimilatory role in A. 



fulgidus (Table Sid). Up-regulated genes; alcohol dehydro- 
genase (AF0024, AF0339) and acyl-CoA transferase/carnitine 
dehydratase (AF0973-AF0974), potentially play a role in co- 
assimilation of organic substrates during growth with T-H2/CO2 
(Table S3). This may also relate to the specific induction of sev- 
eral of the genes of lipid metabolism corresponding to fatty acid 
activation and beta-oxidation (Table Sla, Figure?); acyl-CoA 
synthetase/ AMP- acid ligase II (fadD/ alkK; 9/16 induced homol- 
ogous genes in genome), acyl-CoA dehydrogenase {acd; 4/14), 
enoyl-CoA hydratase (fad; 3/5), 3-hydroxyacyl-CoA dehydroge- 
nase (hdb; 3/10) and acetyl-CoA acetyltransferase (acaB; 3/15). 
The genes of acyl-coenzyme A synthetase/ AMP- (fatty) acid ligase 
(acs; 5/8) may have a role in fatty acid activation or acetate fix- 
ation together with a 3-hydroxy-3-methylglutarate CoA synthase 
homolog (AF0033). The induced genes of putative sterol carrier 
protein (2/3; AFl 174, AF1678) and a short chain fatty acid trans- 
porter (AF1538) may play a role in transport of lipids across the 
membrane. The genes linked to fatty acid oxidation may form 
part of the 3-hydroxypropionate/4-hydroxybutyrate pathway of 
CO2 assimilation (Figure?) (Berg et al, 2007). In relation to 
fixation of CO2, the gene of the large-subunit of ribulose bispho- 
sphate carboxylase (RuBisCo; AFl 587, Table S3) was also highly 
up-regulated during growth with T-H2/CO2. Carboxylase activity 
has been verified in A. fulgidus (Watson et al, 1999). However, the 
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Table 1 | Differential expression (fold change; above 1.5 in bold, ANOVA; p-value 0.00001) of selected genes corresponding to processes in 
Figure 5, graphs display transcriptional abundance (1-3 fold average expression) and standard deviance around mean (dark area). 

■D 

Locus Locus NCBI re O H2 vs. Thiosulfate Late-log vs. 

^ O 

tag annotation (/) o Lactate vs. Sulfate mid-log 

SULFATE REDUCTION fllH^^^^^^^^^^^^^^^^^^^^^^^^feHHHIHHI^^^^^^^^H 



Sulfate + ATP + H2O + 2.0 e" ^ Sulfite + AMP + 2 Pi 



<b N*" \^ \^ 



AF1667 sat 


Sulfate adenylyltransferase 


+ 


P 


-1.16 




AF1668 - 


Hypothetical protein 


+ 


S 


-1.51 




AF1669 aprB 


Adenylylsulfate reductase, subunit B 


+ 


C 




^ [ [ J 


API 670 aprA 


Adenylylsulfate reductase, subunit A 


+ 


C 


-1.40 




AF0756 ppxl 


Putative nnanganese-dependent inorganic 


+ 


C 


-2.39 






pyrophosphatase 










Sulfite + 6.0 e" 


^ Sulfide -1- 3 H2O 










AF0423 dsrA 


Sulfite reductase, subunit alpha 


+ 


C 




mil 


AF0424 dsrB 


Sulfite reductase, subunit beta 


+ 


C 




mil 


AF0425 dsrD 


Sulfite reductase, subunit gannnna 


+ 


C 







AF2228 dsvC idsrC) Sulfite reductase, desulfoviridin-type subunit 
gannnna 



Putative thiosulfate reductases (Thiosulfate + 2.0 e~ Sulfite + Sulfide) 



AF0157 - 


Molybdopterin oxidoreductase, iron-sulfur 


+ 


C 










binding subunit 












AF0158 - 


Hypothetical protein AF0158 


+ 


S 








AF0159 - 


Molybdopterin oxidoreductase. 


+ 


C 










nnolybdopterin-binding subunit, putative 












AF0174 - 


Molybdopterin oxidoreductase, nnennbrane 




C 




-1.23 






subunit 












AF0175 - 


Molybdopterin oxidoreductase, iron-sulfur 




C 




-1.19 


■ ■■■■^ 




binding subunit 












AF0176 - 


Molybdopterin oxidoreductase. 




C 




-1.20 






nnolybdopterin-binding subunit 












AF2384 - 


Molybdopterin oxidoreductase. 


+ 


c 


1.28 








nnolybdopterin-binding subunit 












AF2385 - 


Molybdopterin oxidoreductase, iron-sulfur 


+ 


c 


2.08 








binding subunit 












AF2386 - 


Molybdopterin oxidoreductase, nnennbrane 


+ 


P 


1.95 




— 




subunit 














AF0086 nrtB-1 


Nitrate ABC transporter, pernnease protein 




p 








AF0087 nrtC-1 


Nitrate ABC transporter, ATP-binding protein 




p 








AF0088 - 


Thiannine biosynthesis protein, putative 




p 









AF0092 cysA Sulfate ABC transporter, ATP-binding protein - G -1.21 

AF0093 Cyst Sulfate ABC transporter, pernnease protein - 0 

AF0094 - Molybdate ABC transporter periplasnnic - P 

substrate-binding protein 



AF0638 nrtC-2 Nitrate ABC transporter, ATP-binding protein - P ^: 

AF0639 nrtB-2 Nitrate ABC transporter, pernnease protein - P 

(Continued) 
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Table 1 | Continued 

■D 

Locus Locus NCBI i O H2 vs. Thiosulfate Late-log vs. 

^ O 

tag annotation co o Lactate vs. Sulfate mid-log 

AF0640 - Nitrate ABC transporter, ATP-binding protein, - P -1.50 

putative 



AF1136 - ABC-type nnultidrug transport systenn, ATPase + V -2.09 -1.75 

connponent 

AF1137 - ABC-type transport systenn involved in + R -2.10 -1.77 

nnulti-copper enzynne nnaturation, pernnease 
connponent 

AP1 138 - Predicted nnennbrane protein + S -2.46 -2.01 

AF1561 tauE Predicted pernnease + R HUB 

AF1562 - Uncharacterized conserved protein + S -1.17 : B B B M Wt 

(C) MEMBRANE-BOUND ELECTRON TRANSFER TO TERMINAL ACCEPTOR (BY hdr HOMOLOGS) 

AF0661 qmoO Heterodisulfide reductase, subunit E, putative - C -1.09 BB B BBl 

AF0662 qrnoB' Heterodisulfide reductase, subunit - C BB B BB- 

A/nnethylviologen-reducing hydrogenase, 

subunit delta 

AF0663 qmoA Heterodisulfide reductase, subunit A - C -1.11 BBBBB- 

(hdrA-l) 




dsrMKJOP 



AF0499 


dsra 


Molybdopterin oxidoreductase, iron-sulfur 
binding subunit 


+ 


C 




AF0500 


dsrP 


Molybdopterin oxidoreductase, nnennbrane 
subunit 


+ 


C 


AUBft 


AF0501 


dsrM' 


Nitrate reductase, gannnna subunit, putative 


+ 


C 


ABUft 


AF0502 


dsrK^ 


Heterodisulfide reductase, subunit D, putative 
(Fe-S oxidoreductase*) 


+ 


C 




AF0503 


dsrJ 


Hypothetical protein 


+ 


S 


■ ■■■■ 



dsrMK/hdrDE 



AF0543 dsrK^ Fe-S oxidoreductase* 

AF0544 dsrK^ Fe-S oxidoreductase* 

AF0545 dsrM^ Nitrate reductase gannnna subunit* 



- C 

- c 

- c ABUft 



AF0546 dsrM^ 
AF0547 c/sr/C' 
AF0754 - 
AF0755 hdrDE 

(dsrMn 
(D) HYDROGENASES 



AF1371 
AF1372 



vhtD-1 
vhuA 



AF1373 vhuG 

AF1374 vhuD 

AF1375 hdrB 

AF1376 hdrC 

AF1377 hdrA-2 



Nitrate reductase, gannnna subunit (narl) + C 

Fe-S oxidoreductase* + C 

Hypothetical protein AF0754 - X 2.97 

Heterodisulfide reductase, subunits E and D, - C 4.69 
putative 



Ni,Fe-hydrogenase nnaturation factor* - C 

Methylviologen-reducing hydrogenase, subunit - C 2.61 
alpha 

Methylviologen-reducing hydrogenase, subunit - C 3.55 
gannnna 

Methylviologen-reducing hydrogenase, subunit - C 3.45 
delta 

Heterodisulfide reductase, subunit B - C 4.68 

Heterodisulfide reductase, subunit C - C 6.30 

Heterodisulfide reductase, subunit A - C 8.15 



AF1378 vhtD-2 Ni,Fe-hydrogenase nnaturation factor* - C 23.99 " ■ ■ 

(Continued) 
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Table 1 | Continued 

Locus Locus NCBI w O H2 vs. Thiosulfate Late-log vs. 

^ O 

tag annotation co o Lactate vs. Sulfate mid-log 



AF1379 hydC Quinone-reactive Ni/Fe-hydrogenase B-type 

cytochrome subunit 
API 380 vhtA F42o-non-reducing hydrogenase 

AF1381 vhtG F42o-non-reducing hydrogenase 

(E) LACTATE DEHYDROGENASES (LACTATE ^ PYRUVATE + 2 e-Y 
AF0394 did d-lactate dehydrogenase, cytochronne-type 

AF0395 noxA-2 NADH oxidase 



12.80 



18.66 
53.76 



1.39 
1.82 



AF0506 Fe-S oxidoreductase - C -2.00 

AF0507 FAD/FMN-containing dehydrogenase fused to - C -1.69 

Heterodisulfide reductase, subunit B 



AF0806 


IctP 


l-lactate pernnease 


+ 


C 


-5.81 




AF0807 


lldD 


l-lactate dehydrogenase, cytochronne-type 




H 


-4.57 




AF0808 


did 


Glycolate oxidase subunit 


+ 


C 


-9.04 




AF0809 


lldB 


Heterodisulfie reductase, subunit D, putative 


+ 


C -5.69 


-4.07 




AF0810 


lldO 


Hypothetical protein 


+ 


X 


-5.86 


-M^ : 


AF0811 


lldP 


Hypothetical protein 


+ 


c 


-5.75 




AF0812 




Hypothetical protein 


+ 


X -2.58 


-2.07 





AF0866 glpK 

AF0867 (//c/D') 

AF0868 - 

AF0869 - 
(F) ACETYL-COA 
ACS/CODH 

ACS 

AF0376 cdhE 

AF0377 cdhD 

AF0379 cdhC 
CODH 

AF1100 cdhA-1 

AF1101 cdhB-1 



Glycerol kinase 
Reductase, putative 

Alkyldihydroxyacetonephosphate synthase 
Hypothetical protein 



PATHWAY 



Acetyl-CoA + 5,6,78-Tetrahydronnethanopterin + H2O - 
tetrahydronnethanopterin + CO2 + CoA + Fdred 
Acetyl-CoA + H+ + 5,6,78-Tetrahydronnethanopterin ^ 
Acetyl-CoA decarbonylase/synthase connplex 
subunit gannnna 

Acetyl-CoA decarbonylase/synthase connplex 
subunit delta 

Acetyl-CoA decarbonylase/synthase connplex 
subunit beta 
CO + Fdo. 



C 
C 
C 
S 

-Fdo. 



1.34 
1.34 
1.36 



5-Methyl-5,6,78- 



5-Methyl-5,6,78-tetrahydronnethanopterin + 
C 1.06 




■ ■111 

■nil 



. + 2 H2O ^ CO2 + Fdred WM 

Acetyl-CoA decarbonylase/synthase connplex 
subunit alpha 

Acetyl-CoA decarbonylase/synthase connplex 
subunit epsilon 



1.57 



2.12 



AF2397 cdhA-2 


Acetyl-CoA decarbonylase/synthase connplex 


+ c 


J J J J J 




subunit alpha 






AF2398 cdhB-2 


Acetyl-CoA decarbonylase/synthase connplex 


+ c 






subunit epsilon 







* Annotation from arCOG; (Wolf et al., 2012). 
'Inferred locus annotation. 



role of RuBisCo in anaerobic Archaea may relate to ribulose 1,5- 
bisphosphate recycling, or AMP metabolism, rather than a carbon 
assimilation mechanism (Sato et al., 2007; Estelmann et al., 201 1). 

Genes related to motility (N) were of flagellin (AF1054, 
AF1055, Table S3) and archaeal flagellar biosynthesis (AF0338) 
were induced during growth with T-H2/CO2, indicating a taxis 



response to substrate. The induced genes of signal transduc- 
tion pathways (T) involved genes encoding proteins with poten- 
tial Per-Arnt-Sim (PAS) domains (AF0277, AF0448, AF1045, 
AF1472, AF2420). These may correspond to cellular redox sen- 
sors that have been linked to chemotaxis in Euryarchaeota, but 
also hydrogenase expression in Bacteria (Taylor and Zhulin, 1999; 
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GpIC 

Conserved domains (pfam) FAD binding FAD-oxidase Fer4_8 CCG 



GpIC 



DUF162 



Desulfovibrio vulgaris 
Hildenborough 



I® 0 (5 > 



3027 3028 



:eo> 



Shewanella oneidensis 
MR 1 



Archaeoglobus fulgidus <^ 



AF0805 AF0806 



B U > 



III l!l Ml^*^ 

1783 1782 ^ 1781 

I^^^^HHj/^f other lactate dehydrogenase homologues in A. fulgidus 

1520 1519 1518 



AF0395 AF0396 



AF0807 AF0808 AF0809 AF0811 AF0812 
«AF0808b» AF0810 



Archaeoglobus sulfaticallidus 



AF0868 AF0867 
AF0507 AF0506 



Asulf_00635 00636 00637 00638 00639 00640 00641 nnR49 00643 



FIGURE 3 I Synteny and domain composition between lactate 
dehydrogenases and associated genes. The genes of A. fulgidus 
correspond to a lactate dehydrogenase (lldEFG), D-lactate dehydrogenase 
(did/dld-ll), and lactate pernnease (IctP) which are found co-located in a wide 



range of bacterial species. Genes with sufficient honnology between bacterial 
and archaeal sequences are colored correspondingly. Although, honnology is 
generally low between genes, the donnain connposition for the gene clusters 
is conserved. 



Lenz et al, 2002; Shaw et al, 2009). Other up-regulated genes 
corresponding to intracellular signaling were: histidine kinases 
(AF0893, AF1483), c-AMP binding (AF0971) and universal stress 
protein (AF1526). In addition, genes encoding several putative 
permeases and transporters were up-regulated during growth on 
hydrogen; permeases (AF0121m, AF0123-AF0124), and putative 
proline permease/sodium: solute symporters (AF0965-AF0966, 
AF0969, AF0981-AF0982). 

Fewer genes were highly induced (97 genes, >1.5 fold) 
by utilization of lactate (S-L, T-L) in comparison to uti- 
lization of T-H2/CO2 (Figures IB, 2). These genes were 
functionally enriched in COG's corresponding to metabolic 
processes such as carbohydrate metabolism (G); where car- 
bohydrate kinase (AF1751) and phosphoglycerate mutase 
(AF1752) were induced above 1.5 fold. During growth 
with lactate (S-L, T-L) only a gene with low homology to 
D-lactate dehydrogenase was specifically induced (AF0507, 
AF0506; Table IE, Figure 3). In the inorganic ion trans- 
port (P) category, up-regulated genes encoding ABC-type 
multidrug transporter (AF1136-AF1140, Table IB), and a 
phosphate ABC transporter (AF1356-AF1360), could play 
a role in substrate uptake or sulfate transport. A region of 
unknown, short, DUF2589 related genes (AF0414-AF0417) 
were also specifically up -regulated during growth with 
lactate. 

The genes encoding enzymes of the dissimilatory sulfate 
reduction pathway are affiliated with two COG categories (P 
and C, Table lA). The genes related to reduction of sulfate, i.e., 
APS formation and reduction; sat and aprAB (AF1667-AF1670, 
Table lA) and the membrane-bound QmoABC complex 
(AF0661-AF0663, Table IC) were highly expressed at all times. 



but significantly down- regulated by less than 1.5 fold in relation 
to a shift in energy donor from lactate to hydrogen. The 
inorganic pyrophosphatase (ppx, AF0756, Table lA) which 
is perceived to drive the formation of APS (Peck, 1962) was 
more than 2 fold down- regulated in relation to growth with 
T-H2/CO2. Differential regulation of these genes may be related 
to energy donor (lactate vs. H2) and not terminal electron 
acceptor, as no regulation was observed between S-L and T-L 
samples. 

DIFFERENTIALLY EXPRESSED GENES WITH THIGSULFATE vs. SULFATE 

Transcriptional up -regulation during growth with thiosulfate 
(T-L, T-H2/CO2, Figure IB) corresponded to significant enrich- 
ment of genes in the energy production category (C, Figure 2). 
Genes regulated more than 1.5 fold during growth with 
thiosulfate belong to an operon of putative membrane inte- 
grated periplasmic thiosulfate reductase (AF2384-AF2386). 
This cluster is one of three gene clusters (AF0157-AF0160, 
AF0173-AF0176, AF2384-AF2386, Table 1 A, Figure 4) of 
which protein expression is induced while A. fulgidus utilizes 
(per)chlorate as an electron acceptor (Liebensteiner et al., 
2013). A previous study has linked the gene products of 
AF0157-AF0160 to the twin arginine translocation (Tat) pathway 
(Coulthurst et al., 2012). The molybdenum-binding sub unit, 
encoded by AF2384, contains a similar Tat signal peptide 
(BLASTp, TatP 1.0: http://www.cbs.dtu.dk/services/TatP/, 
Figure 4), indicating a periplasmic location of the thiosulfate 
reductase. 

The second copy of the carbon monoxide dehydrogenase 
gene cdhAB-1; (AFllOO-AFllOl, Table IF) (Dai et al, 1998) 
was unexpectedly highly up -regulated during growth with 
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[4Fe4S] binding 

Conserved domains (pfam) Nitrate_red_gama FAD binding_2 Fer4» CCG FpID 



Vht and Mvh:Hdr 



HdrA 



mmmmm 

\ 1381 1380 1379 1378 1 1377 1376 1375 1374 1373 1372 1371 



Qmo 



HdrDE 



DsrMK-DsrMKK 



DsrMKJOP 



Molybdopterin 
oxidoreductases 



— — dC^o — > 



i>^<^ ^^■^ 



6^^° t^<^ (f><^ e.^^^ 





2384 2385 2386 



00 "Tat signal peptide 



FIGURE 4 1 Synteny and homology between hydrogenases and stripes), Heterodisulfide reductase (CCG donnain; orange), 

menaquinol oxidases potentially involved in energy conservation molybdopterin-binding oxidoreductases; putative 

in A. fulgidus. Hypothetical bifurcating and ferredoxin interacting thiosulfate/polysulfide/thetrathionate or fornnate (AF1202-AF1203) 

genes (green), Menaquinone reductase/Menaquinol oxidase (pink/purple dehydrogenase connplexes (peach). 



thiosulfate. Finally, genes of cobalamin/vitamin B12 biosynthe- 
sis (AF0724-AF0727, AF1843) were specifically up-regulated in 
cultures utilizing thiosulfate. Vitamin B12 is a cofactor in methyl 
transferases, such as the second subunit of the acetyl-CoA trans- 
ferase (ACS)/CODH complex (Banerjee and Ragsdale, 2003). 

Significantly enriched genes up -regulated in the presence of 
sulfate (S-L, Figure IB) belong to a specific set of 5 genes cat- 
egorized under energy production and metabolism (C) corre- 
sponding to a region of one lactate permease and 5 putative 
lactate dehydrogenase (LDH) genes (Ictp; AF0806, lldD; AF0807, 
did; AF0808 and lldEFG; AF0809-AF081 1, Table IE, Figures 3, 5). 
The genes AF0809-AF0811 may encode an oligomeric LDH 
(lldEFG) based on the presence of conserved domains with bacte- 
rial LdlEFG (Figure 3) in Shewanella oneidensis MR- 1 and Bacillus 
suhtilis (Chai et al, 2009; Pinchuk et al, 2009) despite low over- 
all sequence identity (Figure 3). Within the genus Archaeoglobus, 



only A. fulgidus and A. sulfaticallidus are known to couple lac- 
tate oxidation to dissimilatory sulfate reduction (Steinsbu et al., 
2010). A genetic comparison of Archaeoglobales revealed a cor- 
responding genetic region only in these two species, with a 
conserved upstream gene encoding a putative regulatory helix 
turn helix (HTH) motif, and a putative ORF— "AF0808b" homol- 
ogous to the gene Asulf00640 (Figure 3). The putative ORF, 
"AF0808b" may encode a protein with a ferredoxin-binding 
domain. 

Other genes up -regulated in the presence of sulfate belong 
to the COG category defense (V), corresponding to a region 
CRISPR-associated proteins (AF1861-AF1868, Table S3); the 
cmrl-6 genes of RNA guided RNase in Pyrococcus furiosus (Hale 
et al., 2009). The induction of these genes is probably not related 
to growth rate as T-L cultures had a higher specific growth rate 
(Figure lA, Table S3). 
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H2 2H+ 



FIGURE 5 I An overview of the proposed central cellular metabolism of 
A. fulgidus. Gene abbreviations catalyzing reactions are italized; 
underlined — if constitutively expressed; and bold — if high signal intensity 
(>3). Gene abbreviations are colored corresponding to regulation of different 
substrates; H2 (dark blue), lactate (light blue), thiosulfate (8203^"; dark red), 
and sulfate (S04^~; pink). Cartoons of central enzynnes represent a rough 



outline of structural features, subunits are colored correspondingly: 
green — represents potential ferredoxin interacting domains; 
blue — hydrogenases; magenta — menaquinone/quinol interacting subunits, 
orange; heterodisulfide reductase (with [CCG] domains). Gene names 
corresponding to this figure can be found in Table 1 and Table S2 — for genes 
of the acetyl-CoA pathway. 



DIFFERENTIALLY EXPRESSED GENES RELATED TO GROWTH PHASE 

Up-regulated genes influenced by late log growth (Figure IB) 
were distinctly enriched in genes corresponding to information 
storage and processing. In the categories transcription (K) and 
replication (L) the genes were differentially expressed below 1.5 
fold change (Figure 2, Table S3). Genes for cellular processes 
and other COG categories: cell wall/membrane/envelope (M), 
sporulation protein (AF1778. COG - S), histone Al (AF0337; 
COG - L), and cell division inhibitor (AF0696; COG - N); were 
all up-regulated over 1.5 fold in both late-log S-L and late-log 
T-H2/CO2 cultures. 

Genes corresponding to mid-log belonged to the COG cat- 
egories of nucleotide transport and metabolism (F), where 
thymidine phosphorylase (AF1341, AF1342) was highly induced. 
Also, transcription (K) -related genes were significantly up- 
regulated, but at minor fold (<1.5). It was therefore evi- 
dence of distinct functions related to a shift in log phase, 
corresponding to transcription, replication and genes for cel- 
lular processes. These genes include processes such as repli- 
cation, histone modification, and metabolism of nucleotides. 
The functions of genes related to a growth phase specific 
response, was distinct from that of a shift in energy or 
electron acceptors, relating to cellular metabolism in general 
(Figure 2). 



ABSOLUTE ABUNDANCE OF TRANSCRIPTS 

The genes corresponding to KEGG pathways, such as ribosome 
(afu030I0), oxidative phosphorylation (afu00190), aminoacyl 
t-RNA synthesis (afu00970), and methane metabolism 
(afu00680), were enriched amongst genes with high signal 
abundance. As were genes of general metabolic pathways 
(afuOl 100) and biosynthesis of secondary metabolites (afuOl 1 10) 
(Figure SIB). 

In the Oxidative phosphorylation pathway (afu00190), genes 
of the F420H2: quionone oxidoreductase (Fqo) (AF1823-AF1833) 
(Kunow et al, 1994; Briiggemann et al, 2000) and the 
archaeal V-type ATPase (AF1158-AFII68) were continuously 
expressed above 3 fold average expression (Figure SIA). Genes 
for dissimilatory sulfate reduction (sat, ppx, aps, and qmo), 
although differentially regulated, were generally expressed at lev- 
els above 3 fold average expression levels. The genes related to 
terminal reduction of sulfite (SOs^", Table lA); dsrAB (AF0423- 
AF0424), membrane-bound heterodisulfide-reductase:quinol- 
oxidase; dsrMKJOP (AF0499-AF0502), and homologous second 
copy of dsrMK(K) (AF0543-AF0545) were constant and highly 
expressed (Table IC, Figure SIA-I). The notable exception was 
that of dsrD (AF0425) and dsrC transcripts (AF2228), which 
were expressed at lower levels (Table lA, Figure SIA) and were 
both induced in late-log S-L samples. These two genes had a 
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corresponding transcriptional pattern (Pearsons r-value; 0.82) 
despite a large genomic distance (~0.6 Mb). 

Genes relating to all steps of the acetyl- Co A path- 
way were expressed at equivalent and constitutive 
levels (Table S2): CODH/ACS complex (cdhAB-2; 
AF2397-AF2398/c^i/zCD£;AF0376-AF0379), along with N5,N10- 
methylenetetrahydromethanopterin reductase {mtr-1, AF1066); 
the F420 -dependent methylenetetrahydromethanopterin dehy- 
drogenase {mtd, AF0714); methenyltetrahydromethanopterin 
(mc/z; AF1935); formylmethanofuran-tetrahydromethanopterin 
formyltransferase (/tr-i; AF2073 and ftr-2\ AF2207) and the 
formylmethanofuran dehydrogenase {fwdABCD; AF1928- 
AF1935). There was therefore no indication of differential 
regulation relating to the acetyl-CoA pathway of A. fulgidus 
between oxidative (lactate) and reductive (H2/CO2) processes. 
Other genes of central metabolic processes related to acetyl- Co A 
were the constitutively highly expressed pyruvate ferredoxin 
oxidoreductases (porABDG; AF1669-AF1702, Table S3). This is 
in line with a bidirectional role of ferredoxin in oxidative and 
reductive carbon metabolism (Figures 5, 7). 

Constitutively highly expressed genes of electron transport 
flavoproteins {etfA; AF0287 and etfB; AF0286, Table Sib), may 
facilitate reactions relating to the homologous ferredoxin reduc- 
tion by bifurcating butyryl-CoA dehydrogenases/Etf complex in 
Clostridium spp. (Li et al., 2008; Buckel and Thauer, 2013). 

The genes AF2378-AF2380 were constitutively highly 
expressed (Table S3). These genes have previously been linked to 
syntrophic growth in D. vulgaris Hildenborough (Scholten et al., 
2007; Fievet et al., 2011). The proteins encoded by these genes 
have conserved domains NifK (AF2378) and CbiA superfamilies 
(AF2380) (NCBI-CDD) (Rubio and Ludden, 2008), indicating 
a role in the biosynthesis of iron-molybdenum cofactors. These 
may function in the synthesis of molybdopterin oxidoreductases 
that are abundant in the genome of A. fulgidus (Klenk et al, 
1997). Other highly expressed genes of hypothetical proteins are 
AF1617-AF1619 (Table S3), which contain multiple transmem- 
brane helixes. The adjacent associated PAS domain encoding 
gene (AF1620) may indicate an unknown regulatory complex. 

Finally, the genes of rubrerythrin and desulfoferrodoxin (rri, 
rr2, dfx; AF0831-AF0833) were expressed continuously at high 
abundance (Table S3). These are putative oxidoreductases, and 
probably function in elimination of superoxides (Rodrigues et al, 
2005). Genes of cytochrome be oxidase (cydA; AF2296, AF2297, 
Table S3) were highly expressed during all conditions. As in 
Desulfovihrio spp. (Ramel et al, 2013), this may allow cou- 
pling the menaquinone pool to oxygen reduction in A. fulgidus 
(Figure 5). 

DISCUSSION 

In the present work, a model of the energy metabolism in A. 
fulgidus for the utilization of lactate and hydrogen with thiosul- 
fate or sulfate as terminal electron acceptors is presented based on 
transcriptome profiling. 

LACTATE METABOLISM 

Lactate is the "classical" substrate of sulfate reducers, and its 
link to energy conservation in Desulfovihrio has been the subject 



of intense study (Keller and Wall, 2011). Several transcriptional 
shifts were observed in A. fulgidus, involving expression of LDH 
and putative LDH genes (Table IE, Figure 3). Our results indicate 
that during growth with T-L, activity of multiple LDH isozymes 
(Figures 3, 5) may occur in A. fulgidus, as suggested in D. vulgaris 
(Keller and Wall, 201 1 ). When sulfate is used as an electron accep- 
tor, oligomeric LdlEFG may operate together with monomeric 
lldD and did in the oxidation of lactate (Figure 5). The conserved 
"modular" domain composition of the proteins encoded by the 
genes did, the ORF "AF0808" and AF0809, may facilitate a multi- 
meric complex that functions as monomeric homologs encoded 
in other species (Dvu3071, Figures 3, 5). 

The presence of a gene cluster with identical arrangement 
in the lactate utilizing A. sulfaticallidus and A. fulgidus, sup- 
ports a potential role of the LdlEFG in linking lactate oxidation 
with sulfate reduction in A. fulgidus. Acquiring the lldEFG gene 
cluster may have been essential for Archaeoglobales in order 
to perform dissimilatory sulfate reduction with lactate as an 
energy source, potentially via the QmoABC complex to APS 
reductase. 

Oligomeric lldEFG is widely distributed in Bacteria, including 
sulfate-veducmg Deltaproteohacteria (Pinchuk et al, 2009; Pereira 
et al, 2011), but has previously not been identified in Archaea. 
Various functions have, however, been suggested for LdlEFG in 
Bacteria. In S. oneidensis MR-1, the LdlEFG is found to stim- 
ulate the activity of Did- II (Figure 3), indicating a functional 
relationship (Pinchuk et al, 2009). Interestingly, in D. alaskensis 
the LdlEFG is required in syntrophic growth with Methanococcus 
(Meyer et al, 2013). In this model an LdhAB-1 (GplCD) catalyses 
the primary oxidation of lactate, and transfers electrons, possibly 
through thiol/disulfide, to an LdlEFG homologous complex. The 
LdlEFG may transfer electrons to the QmoABC complex, which 
facilitates menaquinol reduction (Meyer et al., 2013). However, 
the LdlEFG is also present in species without a QmoABC complex 
and functions independently as a membrane associated L-LDH 
capable of reducing quinone (Chai et al., 2009; Pinchuk et al., 
2009; Thomas et al., 2011). In order to verify the specific role 
of the LdlEFG homologs in A. fulgidus, biochemical studies are 
required (enzyme activity and protein-protein interaction) to 
understand its relation to Qmo and energy conservation. Perhaps, 
prior to construction of deletion mutants as a genetic system is 
not yet available for this species. 

With the exception of the cdhAB-1 (see next section), genes 
encoding the acetyl- Co A pathway were constitutively expressed 
at high levels (Table S2). This was also true for the F420H2: 
quinone oxido reductase (Fqo) complex, which probably catalyzes 
proton translocation utilizing F420H2 generated by the oxidative 
acetyl-CoA pathway (Briiggemann et al., 2000). The hydroge- 
nases in A. fulgidus were specifically induced during growth with 
hydrogen, and low transcriptional expression of hydrogenases was 
observed during growth on lactate (Table 1). Therefore, it maybe 
questioned whether "hydrogen cycling" (Odom and Peck, 1981; 
Kulkarni et al., 2009) is used as a mechanism for energy conser- 
vation with lactate as the energy source. This would emphasize 
the role of the respiratory Fqo complex and a menaquinone- 
based respiratory system (Figure 5) in energy conservation in A. 
fulgidus during growth with lactate. 
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Several distinct putative menaquinol oxidase:Hdr complexes 
are present in the genome of A. fulgidus. The DsrMKJOP 
(AF0499-AF503) complex and the DsrMK(K) (AF0543-AF0544) 
were constitutively highly expressed (Table IC), whereas a second 
DsrMK (AF0543-AF0544) was expressed at average expression 
levels (Table IC). Multiple membrane-bound DsrMK complexes 
may therefore oxidize the menaquinol (MQH2) generated by the 
Fqo complex (Figure 5). 

The DsrK components may transfer electrons to DsrC by 
breaking the disulfide bonds between the two C — terminal cys- 
teines of the enzyme (Mander et al., 2005; Oliveira et al, 2008; 
Grein et al, 2010). The dsrC gene (AF2228) is, however, expressed 
at average transcriptional abundance vs. high transcriptional 
abundance for dsrAB (Table IC, Figure SI). This is lower than 
previously estimated in D. vulgaris, where the gene of dsrC is 
expressed at high levels (Wall et al, 2008). Although these values 
are more rigorously estimated in the previous study, our results 
point toward a lower expression ratio between dsrAB and dsrC in 
A. fulgidus. This may indicate involvement of additional electron 
transport components from Hdr to DsrAB. However, other elec- 
tron carrying proteins, such as ferredoxin, are expressed at similar 
levels as dsrC (fdx, <1.4 average expression. Table S3). The true 
significance of the role of electron flow via DsrC requires further 
evaluation on translational level. 

REDUCTION OFTHIOSULFATE 

The mechanism of thiosulfate reduction and the import of sul- 
fate for cytoplasmatic reduction is uncharacterized in A. fulgidus. 
The specific growth rate of A. fulgidus cultivated with lactate 
was increased by the utilization of thiosulfate, vs. sulfate, as 
terminal electron acceptor (Figure lA). The reduction of thio- 
sulfate is thermodynamically favorable (Badziong and Thauer, 
1978). However, utilization of thiosulfate vs. sulfate is reported 
as inhibiting for growth rate in D. vulgaris Hildenborough, and 
has been attributed to the toxicity of increasing intracellular 
concentrations of sulfite (Badziong and Thauer, 1978; Pereira 
et al., 2008). The genes corresponding to thiosulfate reduc- 
tase in A. fulgidus are identified by specific up-regulation of a 
molybdopterin oxido reductase (AF2384-AF2386, Table lA). This 
reductase is active during both lactate and H2 -oxidation, and is 
probably a membrane-integrated complex with a periplasmic fac- 
ing active site (Figure 5). The presence of a Tat signal peptide 
(Figure 4) indicates that export is facilitated by the twin arginine 
translocation pathway (Coulthurst et al, 2012). 

A periplasmic reduction of thiosulfate may exclude the build- 
up of toxic intracellular levels of sulfite and may partly explain 
the high growth rate observed for A. fulgidus during cultivation 
with lactate and thiosulfate. It is unlikely that the reduction of 
thiosulfate to sulfite (E°^ —402 mV) contributes to energy conser- 
vation (Badziong and Thauer, 1978; Stoffels et al, 2012). Rather, 
final intracellular reduction of SOa^" to S^~ (E°^ — 116 mV) has a 
redox potential sufficient for energy conservation (Thauer et al., 
2007). In most Desulfovibrio spp., an indistinguishable ion gradi- 
ent symport has been found for thiosulfate and sulfate (Cypionka, 
1987; Stahlmann et al., 1991). However, such a mechanism has 
not been identified in A. fulgidus (Rabus et al., 2006). The 
genes previously annotated "sulfate ABC transporter permease" 



(AF0092-AF0094) are probably a molybdate-specific transporter 
(Klenk et al, 1997; Hollenstein et al, 2007). The induction of 
genes corresponding to a ABC-type transport system (AF1136- 
AF1138) during growth with sulfate (S-L) may ambiguously be 
assigned as a putative sulfate transporter, as these genes are also 
induced during growth with lactate (Table IB, Figure 5). The 
gene tauE is proposed to encode a sulfite exporter in Cupriavidus 
necator {Ralstonia eutropha) during sulfoacetaldehyde degrada- 
tion (Weinitschke et al., 2007). The constitutive highly expressed 
tauE homolog (AF1562) may be assigned a putative function for 
sulfite import in A. fulgidus (Table IB, Figure 5). The utilization 
of thiosulfate is a common property of all Archaeoglobus spp. 
and Ferroglobus placidus, however, homologous of the putative 
periplasmic AF2384-AF2386 gene cluster can only be found in the 
species A. fulgidus and F. placidus (BLASTp, Absynte, and Syntax- 
tools Oberto, 2013). The DsrAB of A. fulgidus displays a high level 
of in vitro thiosulfate reductase activity (Parey et al, 2010) and 
may play a role as a parallel process of cytoplasmatic thiosulfate 
reductase. Although a common trait, different Archaeoglobus spp. 
seem to have diverging enzyme systems for thiosulfate reduction. 

Unexpectedly, when thiosulfate was substituted for sulfate as 
electron- acceptor, a second copy of the cdhAB-1 subunits in the 
ACS/CODH complex was induced (AFllOO-AFllOl, Table IF, 
Figure 5). The specific regulation corresponding to terminal elec- 
tron acceptor (thiosulfate) may indicate a preferential utiliza- 
tion of different electron carriers between dissimilatory sulfate 
and thiosulfate reduction. A similar shift in genes of cobal- 
amin/vitamin B12 biosynthesis may also affect the function of the 
ACS sub unit (Banerjee and Ragsdale, 2003). Previous studies have 
shown that regulation of CODH/ACS complexes in M. acetivorans 
are related to carbon source (Matschiavelli et al., 2012). 

HYDROGEN METABOLISM 

Archaeoglobus fulgidus possesses only two hydrogenases; the 
periplasmic Vht hydrogenase and the soluble Mvh:Hdl (Mander 
et al, 2004). The latter may take part both in energy conservation 
and in generation of Fdred for CO2 — fixation through the acetyl- 
CoA pathway. The reductive acetyl- Co A pathway requires at least 
3 mol Fdred for the generation of one mole pyruvate from CO2 
(Fuchs, 2011). Similar to the methanogens, a bifurcation reac- 
tion is obligate in A. fulgidus for the generation of Fdred while 
growing autotrophically with hydrogen. In addition, Fdredcan be 
hypothesized as an electron donor to APS reductase, through the 
QmoABC (Ramos et al, 2012) and the DsrAB (Oliveira et al, 
2008, 201 1) offering a potential coupling between ferredoxin and 
electron transport phosphorylation. 

The genes of periplasmic Vht hydrogenase represented the 
highest transcriptional shift of any genes in relation to growth 
on T-H2/CO2 and was expressed at a high level relative to 
average signal abundance (Table ID, Figures 4, 5). The result- 
ing two protons from a periplasmic hydrogenase reaction cat- 
alyzed by Vht may contribute directly to the formation of pmf 
during growth. In Methanomicrobiales, the Vht hydrogenase 
homolog donates electrons via methanophenazine (MP) to a 
cytoplasmic- facing, membrane-bound Hdr (HdrDE, Figure 6A) 
(Deppenmeier et al., 1992; Ide et al., 1999; Thauer et al., 
2010). Similarly, the VhtABC complex in A. fulgidus may reduce 
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menaquinone (MQ). A subsequent menaquinol (MQH2) oxi- 
dation, facing the periplasm, by a membrane-bound Hdr may 
translocate two protons (Figures 5, 6A). The observed co- 
induction of vht hydrogenase genes and a fused hdrDE homolog 
encoding dual [CCG] domains (AF0755, Figure 4) suggests a 
close physical interaction between the two encoded complexes 
that may form a distinct path of electron flow to DsrAB. However, 
the genes of DsrMKJOP or DsrMK(K) were constitutively 
expressed, and therefore, electron flow is also possible via these 
complexes (Table 1, Figures 4, 5). The reaction probably repre- 
sents the major pathway of energy conservation during growth 
with H2. 

In contrast to most methanogens and SRP, A. fulgidus pos- 
sesses only one potential mechanism for ferredoxin generation 
from hydrogen; namely the Mvh:Hdl(Hdr) catalyzed reaction 
(Klenk et al, 1997; Thauer et al, 2010; Pereira et al, 2011) 
(Figure 6A), The Mvh:Hdr hydrogenase is one of, so far, 4 per- 
ceived reaction mechanisms for the reduction of ferredoxin from 
H2 during autotrophic growth (Fuchs, 2011) (Figures 6A,B). In 
methanogens the Mvh:Hdr is the key enzyme of the recently 
named Wolfe cycle (Thauer, 2012), which catalyses the cru- 
cial bifurcation reaction that couples the first (Fdred is required 
for the fixation of CO2) and last step of methanogenesis 
(reduction of heterodisulfide, CoM-S-S-CoB). No net Fdred is 
generated from this reaction, as generation of Fdred ^nd het- 
erodisulfide reduction are interdependent (Figure 6A). In order 
to assimilate carbon — an anaplerotic hydrogenase — the energy- 
conserving membrane-associated hydrogenase (Ech) is required 
in these methanogens for additional generation of Fdred for 
anabolic processes (Figure 6B) (Lie et al, 2012; Thauer, 2012). 

The presence of ferredoxin-binding sites ([4Fe-4S] clusters) in 
the structures of DsrAB indicate that soluble Fdred or a ferre- 
doxin reductase complex may facilitate the steps of two -electron 
transfer to S03^~; from redox state +IV, to +11 and 0 (Schiffer 
et al, 2008; Oliveira et al, 2011). However, if the reduction of 
sulfite required 2 mol Fdred prior to reduction by I mol reduced 
DsrC, the pool of available oxidized DsrC would soon be depleted 
(Figure 6A). Therefore, disulfide (X-S-S-X) would not be avail- 
able for recycling Fdox to Fdred by Mvh:Hdl mediated bifurca- 
tion.If the reduction of sulfite was dependent on only 1 Fdred? 
anaplerotic hydrogenase would still be required for the genera- 
tion of Fdred for subsequent CO2 fixation; analogous to the Wolfe 
cycle (Thauer, 2012). During growth with hydrogen, the absence 
of an anaplerotic ferredoxin reductase in A. fulgidus requires mul- 
tiple two-electron transfers for the reduction of sulfite by other 
mechanisms, either by an unknown electron donor or repetitive 
association, oxidation and dissociation of DsrC. 

According to this model Fdred is not a viable electron donor 
for reduction of sulfite in A. fulgidus during growth with H2. 
The electrons for reduction of sulfite must therefore be provided 
by the Vht hydrogenase (Figure 6A). Hence, Fdred is probably 
utilized in biosynthesis rather than energy conservation during 
growth with T-H2, and may be a plausible explanation to the low 
transcriptional levels of mvhihdl (Table ID, Figure 6A). 

Requirement of a Fdred -driven "confurcation" mechanism via 
the QmoABC complex for APS reduction in A. fulgidus would 
according to our model inhibit fixation of CO2 (Figure 6). 



Accordingly, no growth with sulfate and hydrogen (S-H2/CO2) 
has been observed for A. fulgidus (Stetter et al., 1987; Steinsbu 
et al, 2010). While growth on sulfate with hydrogen (S-H2) 
does not occur; A. fulgidus is capable of utilizing sulfate as 
terminal electron acceptor with CO or formate as electron donors 
(Stetter et al, 1987; Henstra et al, 2007). The redox potential 
of CO/CO2 indicates the capacity to reduce ferredoxin directly 
without the need for a bifurcation reaction (Thauer et al, 2007). 
The redox potential of formate/C02 is similar to that of hydrogen 
(Thauer et al., 2007), and requires bifurcation for the genera- 
tion of Fdred- In addition, A. sulfaticallidus grows on S-H2/CO2 
(Steinsbu et al, 2010) and genome analysis of A. sulfaticallidus 
(Stokke et al., 2013), did not provide an alternative mechanism of 
Fdred generation. 

An unknown cytoplasmic formate dehydrogenase (AF1203- 
AF1202, Figure 4; Henstra et al., 2007) probably associates with 
the HdrA subunit of the Mvh/Hdl complex in order to catalyze 
the reduction of ferredoxin (Lie et al, 2012). In order to escape 
the proposed physiological impasse of Fdred as an intermediate of 
APS reduction, a formate dehydrogenase would also be needed 
to associate with QmoB (a HdrA homolog. Figure 4) and drive 
energy conservation by a confurcation mechanism. An unknown 
mechanism may also be facilitated by the gene product of AF 1238 
(Figure 4). Similarly, in A. sulfaticallidus, the Mvh hydrogenase 
subunit may, tentatively, form two complexes; one with HdrA and 
one with a homologous QmoB. 

Considering growth with formate and sulfate, and the sim- 
ilar genomic composition of A. fulgidus and A. sulfaticallidus; 
the most plausible explanation for the inability of A. fulgidus to 
grow on S-H2 is a regulatory link between hydrogen and observed 
down- regulation of a pyrophosphatase gene {ppx, AF0756). 
Additionally, we observed minor down- regulation of genes cor- 
responding to the Sat-ORF2-AprAB operon (Table lA, Figure 5). 
A reduced expression of Ppx would inhibit or limit the formation 
of APS, as the pyrophosphatase reaction drives the total reaction 
to completion (Peck, 1962). Uniquely for A. fulgidus, the ppx 
gene is close to the inversely induced membrane-bound Hdr gene 
{hdrDE; AF0755) located on the opposite strand (Tables 1A,C, 
Figure 4), which may suggest a regulatory link. Therefore, the 
inability of A. fulgidus to grow with sulfate and H2 may 
relate to transcriptional regulation rather than a physiological 
limitation. 

Clearly, further biochemical characterization is needed to ver- 
ify the proposed regulatory mechanism. There is also a need to 
characterize a mechanism for growth with sulfate and formate for 
A. fulgidus, and S-H2 for A. sulfaticallidus (Figure 6A). 

GENERATION OF F420H2 IN THE ABSENCE OF A DEDICATED 
HYDROGENASE (Frh) 

Archaeoglobus fulgidus lacks the cytoplasmic F420 -reducing hydro- 
genase (FrhABG) that catalyzes the reduction of F420 in most 
methanogens (Alex et al, 1990; Thauer et al, 2010) (Figure 6C) 
Therefore, mechanism for generating the reduced F420H2 
required for carbon fixation through the reductive acetyl- Co A is 
unknown (Figure 5). A negligible role of NAD(P)H is supported 
by low expression levels of F42oH2:NADP^ oxido reductase genes 
in our study (AF0892, AF1209; Table S3). 
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FIGURE 6 I (A) A schematic comparison of hydrogenotrophic and 
methanogenic metabolism with that of A. fulgidus. The Wolfe cycle of 
autotrophic hydrogenotrophic nnethanogens without cytochronnes is shaded 
in red. The Mvh:Hdr is the only known ferredoxin-reducing connplex present 
in A. fulgidus. A pathway analogs to the Wolfe cycle is indicated in orange — if 
Fdred is required for either APS or SOs^" reduction, prior to the reduction of 
di-thiol (X-SH HS-X), an anaplerotic ferredoxin-reducing hydrogenase is 
required. (B) An overview of alternative complexes for Fd reduction (dotted 



borders), which are absent in A. fulgidus. Most nnethanogens and SRP 
nnaintain several of these complexes. The Wolfe cycle requires an anaplerotic 
ferredoxin (Fd)-reducing hydrogenase (Ech/Eha), in order to fixate carbon in 
anabolic processes. (C) Putative mechanisms for generation of F420H2. 
Complexes absent in A. fulgidus outlined in (B): the electron bifurcating 
[Fe-Fe] hydrogenase (Huang et al., 2012), the Rnf complex (Biegel et al., 2011 ; 
Tremblay et al., 2013), the energy converting hydrogenase (Ech), (C): the 
F42o-reducing hydrogenase (Frh) (Thauer et al., 2010). 



Independently of the Fpo complex and Frh hydrogenase; the 
FpoF subunit is shown to reduce F420 coupled with oxidation of 
Fdred in M. mazei (Welte and Deppenmeier, 2011). It is therefore 
possible that FqoF in A. fulgidus catalyzes a similar mechanism for 
the generation of reduced F420 (Figure 6C, /^oF: AF1833, Table 
S2). The required Fdred must be provided by the bifurcation reac- 
tion facilitated by the Mvh:Hdl hydrogenase (Figures 6 A,B). As 
discussed in the previous section, the main route of energy con- 
servation is probably provided by the periplasmic hydrogenase. 
The low transcriptional levels of mvhihdl may be a reflection 
of translational levels of Mvh:Hdl hydrogenase, if an alternative 
pathway of F420 reduction is present independently of Fdred- 

Vorholt et al. (1995) suggested the possibility that reduced 
F420H2 may be generated by reverse electron flow through 
menaquinol oxidation. The Fqo complex, including the FqoF 
subunit, is also a potential MQH2 oxidase (Figures 5, 6C). 
The constitutive expression of the entire Fqo complex indi- 
cates that the complex may function in reverse as a pmf 
(|xAH+) -dependent menaquinol — F420 oxido reductase, where 
the menaquinol (E°^ — 75mV) may donate electrons for the 
reduction of F420 (E°^ — 360mV). The resulting positive redox 
potential (E°^ +285 mV) may be overcome in a process assisted 
by the consumption of proton gradient. Further support for 
such a mechanism can be found in the common ancestry of 
the respiratory complex I and Ech hydrogenase ( Redder ich, 
2004; Moparthi and Hagerhall, 201 1). In mitochondria and iron- 
oxidizing Thiobacillus ferrooxidans the NAD(P)H-oxidoreductase 
(Complex I) has been shown to catalyze this reaction at the 
expense of ATP hydrolysis, which is perceived to be coupled 
to generation of a pmf by reversal of ATPase (Chance and 
Hollunger, 1960; Vinogradov, 1998; Elbehti et al, 2000). The 



Eha/Ech hydrogenase activity is linked to the reduction of ferre- 
doxin (Figure 6B) (Meuer et al., 1999), as the reverse electron 
flow of the Ech dehydrogenase catalyzes the formation of CO 
from CO2, and H2by consumption of pmf (Bott et al, 1986; Bott 
and Thauer, 1987; Lie et al., 2012). Therefore, the reduction of 
Fdred(E°' -500 mV) with H2 (E°^ -414 mV; -300 mV at 10 Pa 
H2) results in a positive redox potential (E°0 of at least +86 mV 
(or +200 mV at 10 Pa H2) and is considered possible with the 
utilization of a proton gradient (Figure 6B) (Thauer et al., 2007). 

It remains to be shown in A. fulgidus whether it is possible to 
drive the reduction of F420 by MQH2 oxidation (E°^ +280 mV) 
and a proton gradient by e.g., constructing deletion mutants of 
fqo genes encoding MQ interacting components, or biochemical 
characterization by inverted vesicles (Baumer et al., 2000) coupled 
to ATP hydrolysis. 

CO-ASSIMILATION OF ORGANIC SUBSTRATES 

The ambiguous roles of genes corresponding to fatty acid 
metabolism (COG; I, Figures 2, 7) during growth with H2 is 
discussed. Putative co-assimilation of organic acids is consid- 
ered more likely than a homologous secondary carbon fixation 
pathway, although, both mechanisms may contribute to conserve 
Fdred during autotrophic growth with H2. 

The genes of methylmalonyl metabolism (AF2215-AF2219, 
Table Sib, Figure 7) were continuously highly expressed. These 
enzymes may serve as a pathway of propionate degradation 
(Takaki et al, 2010; Moon et al, 2012). During growth with 
T-H2/CO2, several genes related to fatty acid biosynthesis 
and metabolism were induced (COG; I, Figure 2, Table Sla). 
Expression of the genes may be affected by trace amounts of 
fatty acids present in yeast extract (0.03% weight) and points to a 
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FIGURE 7 I Genes corresponding to fatty acid and propionate 
metabolism, compared with putative steps in the 4-hydroxybutyryl (4HB) 
and 3-hydroxyproprionyl (3HP) pathway of autotrophic Creanarchaea. 

Honnologous genes correspond to several steps in the 4HCD and 3HP pathway 
(internnediate nannes in gray), and nnay not be easily-distinguished by honnology 



alone. Genes of the acetyl-CoA pathway are expressed at a higher level than the 
key enzynnes of putative 4-hydroxybutyryl-CoA dehydratase (4HCD), which are 
not unifornnly induced. Genes that are not induced by T-H2/CO2 are underlined, 
for up-regulated genes — the number of up-regulated vs. total number of 
homologs is indicated (genes are annotated in Table SI ). 



potential for fatty acid scavenging/co- assimilation of organic sub- 
strates during autotrophic growth (Klenk et al, 1997; Zarzycki 
and Fuchs, 2011). 

Recently, Parthasarathy et al. (2013) demonstrated induced 
activity of phenylalanine degradation in the presence of yeast 
extract. Despite amino acids being a major component of yeast 
extract, none of the putative genes reported in the previous study 
(Table 2 in Parthasarathy et al, 2013) were induced during T- 
H2/CO2 growth in this study. Thus, there is no uniform induction 
of putative scavenging mechanisms for organic carbon during 
autotrophic growth. 

The genes related to fatty acid biosynthesis also encode 
enzymes in the 3-hydroxypropionate/4-hydroxybutyrate 
(3HP/4HB) cycle identified in Metallosphaera sedula, and could 
potentially represent a secondary carbon fixation pathway in A. 
fulgidus (Berg et al, 2007; Estelmann et al, 2011). The genome 
of A. fulgidus includes 3 homologs of the 4-hydroxybutyryl-CoA 
dehydratase {4hcd), (AF0333, AF0885 and AF1027, Table Sic), 
which is a key enzyme of the 4-hydroxybutyrate carbon dioxide 
assimilation pathway (Figure?) (Berg et al, 2007). During 



growth with T-H2/CO2, one of the 4hcd homologs (AF0885) 
was induced (1.5 fold). The induced gene displayed average 
signal intensity (1.3), and was expressed at a similar level as the 
two other unregulated homologs (AF0333 and AF1027). The 
differential transcriptional expression of 4hcd (Msedl321) in 
M. sedula was related to autotrophic vs. heterotrophic growth 
and resulted in a more than 7 fold up -regulation (Auernik and 
Kelly, 2010). The presence of a 3HP/4HB cycle was refuted by 
Estelmann et al. (2011) by a subsequent study on the obligate 
autotroph "A. litotrophicusl' where enzyme activity of key 
processes could not be detected. The presence of 5 homologs of 
A. fulgidus 4hcd in the genome of Desulfatihacillum alkenivorans 
suggest that this enzyme is involved in alkene degradation in 
these species (Estelmann et al., 201 1). Analogously to A. fulgidus, 
the facultative autotroph and fatty- acid and alkene degrading 
D. alkenivorans utilizes the bacterial acetyl- Co A/ Wood- Ljungdal 
pathway (So and Young, 1999; Callaghan et al, 2012), indicating 
physiological similarities between the distantly related species. 

Alternative mechanisms for up -regulation of homologous of 
propionate and beta- oxidation may be co- assimilation of organic 
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substrates that may supplement the reductive acetyl-CoA path- 
way. In the hght of the considerations of the role of Fdred during 
growth with T-H2/CO2, this may provide a significant advantage 
by supplementing reduction of CO2 with the uptake of reduced 
organic acids. The constitutively highly expressed EtfAB may also 
provide a source of Fdred by an unknown bifurcation reaction 
from ambient fatty or amino acids (Buckel and Thauer, 2013; 
Parthasarathy et al, 2013). In summary, our data may add sup- 
port to the theory that the 3HB/4HP cycle may have originated 
from a heterotrophic pathway; or as a co-assimilatory pathway in 
Archaea (Fuchs, 2011; Zarzycki and Fuchs, 2011). 
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Figure S1 | Distribution of array abundance and correlation with and 
central metabolic processes. (A) Distribution of intensity of 
quantile-nornnalized data; (A-l) distribution of values corresponding to 
nnininnunn expression per gene and (A-ll) nnaxinnunn expression per gene. 
(B) Gene set enrichnnent analysis (GSEA) (Subramanian et al., 2005) on 
the distribution of selected genes on a list of genes sorted by nnininnunn 
intensity values (i, corresponding to pie chart A-l), i.e., deviation fronn a 
randonn distribution, is displayed for identified peptides and genes 
corresponding to KEGG pathways. The nnaxinnunn enrichnnent score (ES) 
corresponds to the largest deviance fronn randonn distribution, the region 
above the nnaxinnunn enrichnnent score corresponds to the leading edge 
(LE), this region is colored corresponding to the intensity at point of ES. 
The nunnbers denote — [values below LE/above LE]. 
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